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Abstract. - We propose a model of nodal order parameters for the superconducting A and B - 
phase of skutterudite PrOs4Sbi2 and discuss the associated angular dependent magnetothermal 
conductivity at low temperatures. The model for the hybrid gap functions A(k) containing an 
s- and g- wave parts is consistent with recent thermal conductivity experiments in PrOs4Sbi2. 
In particular the model accounts for the data on polar and azimuthal field angle dependence 
of Kzz{9, (f>). The low temperature behaviour of thermodynamic properties in zero field is also 
presented. We show that the effect of impurity scattering on a nodal hybrid gap function 
immediately leads the opening of a gap and related exponential behaviour of low temperature 
specific heat and thermal conductivity which is very different from d- wave superconductors. 



Introduction. - Superconductivity in the cubic (T^ symmetry) Heavy Fermion (HF) 
skutterudite PrOs4Sbi2 with K has provoked great interest since it exhibits a number of 

characteristics that suggest the presence of nodes. |1I2I3| However there is also a proposal that 
it is a conventional s-wave superconductor. {4; The specific heat results |1I2| indicate that firstly 
a low temperature power law behaviour and secondly the presence of two superconducting (sc) 
phases with two consecutive specific heat jumps at rci = l-82 K and Tc2=1.75 K somewhat 
similar to the observations in UPts. More recent low temperature thermal conductivity results 
|S] in the vortex state confirm i) the presence of two superconducting A (T < Td) and B 
(T < Tc2) phases with different nodal structure, ii) the presence of point nodes in both A 
{H > 0.75T, T < Tc) and B {H < 0.75T, T < Tc) phase. More specifically the thermal 
conductivity indicates the presence of four point nodes at k= (1,0,0), (0,1,0), (-1,0,0) and 
(0,-1,0) (in units of kp) in A(k) in the A-phase while in the B-phase the two point nodes are 
located at k= (0,1,0) and (0,-1,0). Present investigation also excludes a fully symmetric A- 
phase order parameter with additional point nodes along [001] direction. 
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Symmetry classification of SC order parameter. ~ The presence of two sc phases is 
possibly connected with the vicinity of a field induced antiferroquadrupolar ordered phase 
above 4T which is likely due to a level crossing of the tetrahedral CEF states P, where one of 
the excited triplet states crosses the singlet ground state Fi. In U-based unconventional 
HF superconductors an exchange of (dipolar) spin fluctuations in the itinerant 5f-electrons is 
frequently implied as the mechanism for unconventional SC pair formation. In the present 
case a new interesting possibility for pair formation is the exchange of inelastic quadrupolar 
Fi- F5 fluctuations of the essentially localized 4f electrons of Pr. Since the AFQ critical field 
is of the same order of magnitude as the sc upper critical field of 11^2 (0) = 2T impeding CEF 
level crossing as function of applied mangetic field will strongly affect the pair potential and 
may lead to a change of the sc gap symmetry. This does however not explain the origin of the 
zero field Tc splitting. Details of this microscopic mechanism and its prefered gap symmetries 
are not clear, hence we have to consider possible order parameters compatible with the above 
node structure on a phenomenological basis. The gap function may be expanded in terms of 
basis functions i/jf (k) which transform like representations F of the crystal symmetry group 
(i=l-d is the degeneracy index, the index 1 denoting the degree of F is suppressed). So far 
there is no information from NMR Knight shift or IIc2- Pauli limiting whether PrOs4Sbi2 
has spin singlet or triplet pairing, although recent /iSR measurements |7] seem to suggest the 
presence of condensate magnetic moments. Only the singlet case will be considered here. The 
gap function should then be given by 

A(k) ^r;fVf(k)^A/(k) (1) 

T,i 

where the form factor /(k) is normalized to one and A is the temperature dependent 
maximum gap value. In the spirit of the Landau theory only a single representation with the 
highest Tc should be realized and for T >Tc, the free energy may then be expanded in terms of 
possible invariants of the order parameter components W Vi which are determined by Landau 
parameters ct^{T) and pf . The node structure is then fixed by the specific symmetry class of 
A(k) defined by the set of rjf . However we are interested in the complementary temperature 
range T -^T^ relevant for magnetothermal conductivity experiments. In this regime Landau 
expansion is unreliable. Especially it is unjustified to assume that only a single representation 
F will have appreciable amplitude in Eq. JQ), one should expect that at low temperatures 
A(k) will be a hybrid gap function, i.e. a superposition of basis functions belonging to a few 
favorable representations which have approximately equal T^- A striking realization of such a 
hybrid gap function has recently been found in the nonmagnetic borocarbide superconductor 
YNi2B2C and possibly also in LuNi2B2C. There A(k) =A/(k) is a sum of two fully 
symmetric (for D4/1) gap functions of different degree ('s+g-wave') given by 

/(k) = \[l-{kt + kl-Qklkl)] (2) 

For hybrid gap functions nodes occur only for special 'fine tuning' of amplitudes, they are 
not enforced by the gap symmetry. ^1^]. Microscopic model calculations ^21 show that 
this is possible for a wide variety of plausible pairing potentials. The case of a nodal hybrid 
order parameter is realized in YNi2B2C to an astonishing accuracy. This is not yet clear 
for PrOs4Sbi2 but the nodal structure refered above cannot be realized with gap functions 
corresponding only to a single Tft,-representation. Experimental evidence from the field angle 
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Fig. 1 Fig. 2 



Fig. 1 - Quasiparticle DOS for the A- phase gap functions (full line) and B- phase gap function 
(dashed line). 

Fig. 2 - Temperature dependence of specific heat Cs (Cn = ^nT) for the A- phase gap functions (full 
line) and B- phase gap function (dashed line). The specific heat jumps ACs/Cn a.t T = Tc are given 
by 0.92 and 1.18, respectively. 



dependence of Kzz{d,4>) leads us to the following simple proposals for hybrid gap functions 
A(k) =A/(k) for PrOs4Sbi2 with four and two nodal points respectively: 



A-phase: /(k) = 1 - A:^ - /c^, B-phase: /(k) = 1 - kfj. (3) 

Both A- phase and B-phase are described by hybrid gap functions which are superpositions 
of three T/i representations. The B-phase gap function has lower symmetry and both are 
threefold degenerate. An alternative gap function for the A-phase of PrOs4Sbi2 would be the 
same model of eq. Q as used already for the borocarbides. 

Quasiparticle DOS and thermodynamics. ~ In the normal state the estimated linear 
specific heat coefficient is Q 7 >350 mJ/mole K^, considerably larger than the value obtained 
from the effective dHvA masses (m* ^2.4-7.6 m). jl3j However in comparison with other Pr- 
compounds both experiments support the HF character of PrOs4Sbi2. We have calculated 
the quasiparticle DOS N(E) for the gap functions of A and B phase which determines the 
low temperature thermodynamics. The DOS functions are shown in Fig. H For E/A <Cl 
the lowest linear term N(E)/No=a(E/A) has a = j for A phase and -I for B phase. Here 
Nq is the normal state DOS. Correspondingly, the low T specific heat, spin susceptibility and 
superfluid density are given by, respectively. 



«^C(3)(?), ^ = a(21n2)^ 



7„r 7r2 A ' Xn ^ 

A-phase: pf(T) pC(^T) = 1 - ca{^)' (4) 

B-phase: p^(T) = 1 - 3^ p^iT) ^ p';{T) = 1 - cb{^Y 

Xn A' 

Here pl{T) = p\{T)/ p\{0) for the axis directions i=a,b,c and ca.b are numerical constants. 
Furthermore the specific heat in the whole temperature region is calculated and given by 
Fig. 13 where we have assumed phenomenologically the pairing interaction V\^\^i to be of the 
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separable form V f {k) f (k' ) . The Cs/i'fnT) vs. T/Tc is found to be nearly independent of V 
for each gap model. However for fixed V the obtained Tc values for A and B phase are quite 
different. We consider now thermodynamics in the vortex state and in the low temperature 
limit with F <C T ^ v^/eH <C A(0) where where v is the (isotropic) Fermi velocity and F the 
quasiparticle scattering rate and vV eH is the typical magnetic energy. The field H is applied 
in an arbitrary direction defined by polar and azimuthal angles 9 and (j) respectively refered to 
a cubic [001] axis. In IHlFlS) it has been established that the dominant effect of the magnetic 
field in nodal superconductors is the Doppler shift in the quasiparticle spectrum introduced 
by the supercurrent around each vortex. A short derivation is also sketched in Ref. 110'. In 
the low T limit the quasiparticle DOS for E=0 in the two phases (i=A,B) is given by 



1\q 4 TT A 

Here a; = v • q/A is the normalized Doppler shift energy and the averaging over the 
Fermi surface and vortex lattice has to be performed. For brevity we use xa = {\x\)a and 
xb = (|a;|)_B. The angular dependent functions 1a. 3(9,4') are given by 

lAi9,(j)) = ^[{1 ~ sin^ 9 sin^ (j))^ + [1 - sin^ 9 cos^ , 
Ib{9A) = i(l-sin2(?sin2 0)^] (6) 
This leads to a field angle dependent specific heat for the A and B phase which is given by 



C,/7nr=-^/.(0,0) . (7) 

Similar angular dependences governed by Ia,b{9, (f>) niay be obtained for the spin suscep- 
tibility Xs(T) and superfluid density PsiT)- In principle the angular dependent specific heat 
is adequate to identify the point nodes in A(k) . Recently a related experiment is reported 
for YNi2B2C. 16 The observed cusps in the angular dependence of Cg clearly indicate the 
presence of point nodes along the tetragonal plane axes at (1,0,0), (0,1,0) etc., which is fully 
consistent with thermal conductivity data [Hj. 

Impurity scattering and magnetothermal conductivity. - In order to compute the thermal 
conductivity, it is imperative to consider the effect of impurity scattering. As already discussed 
elsewhere the impurity scattering in "s+g"- wave supercondutors is completely different from 
other nodal superconductors without s- wave component |17I18) . Since the unitary limit gives 
practically the same result as Born limit, we can follow Abrikosov and Gor'kov and we 
have 



Cul + (A + Af 



l\2 



(8) 
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where ujn and A are renormahzed Matsubara frequency and so order parameter respec- 
tively. Here a = | and f'= |-k!^-k^ for A phase while a = | and f'=i-ky for the B phase. 
First we consider H=:0. In the limit of small F it is shown that A= aA+F which implies that 
a quasiparticle energy gap opens up immediately. These energy gaps are well approximated 

by ujg ~ F(l + ^y^- This leads to C,./7„T and ku/T - (^) ^e"'^''/'^. Therefore unlike 
d- wave superconductors there will be no universal heat conduction |2()ll21j . The thermal 
conductivity vanishes exponentially when T /ujg <ti 1. 

In the presence of a magnetic field ujn in r.h.s of eq. ((HJ have to be replaced by 
iv -q where v -q is the Doppler shift, then for ^ ^1 the quasiparticle damping Co=limi^^o ^ 
in A and B phase is given by 



where x^ (i=A,B) are defined below eq. |SJ| and b^i = 2, bs = 4. A similar expansion in 
terms of x; up to second order is used to obtain the thermal conductivity tensor: 



21n^r2;A^ l + %xl ^ ' 



where we have defined x' = 1 ^'■^^ (1 — sin^ cos^ 2 . The component Kyy may be 
obtained by replacing x' —f xb in n^^. Similarly for the B- phase we obtain 



VV 



XB 



1 + ^2^B + llx'^g 



ln(— ) 

\ Xb / 

3 1 



3tt 



XB 



I 31^2 
5 ^2 



21n(^) 

\xb / 



1 + 31 2 
^ + 64 "^-B 



(11) 



Note that in k-^^ and Kyy have no linear term in Xi and therefore don't exhibit a VH 
behaviour. Indeed eqs. reproduce the (j) dependence of nf^^ observed in JF. For 9 = 

tt/2 cusps in Hzz at (j) — 0, tt/2 (A phase) and = 7r/2 (B phase) are predicted in accordance 
with experiment. On the other hand the present theory predicts the dependence of Kzz 
while the data appear to be approximately linear in H. Also when the heat current is parallel 
to a pair of nodal directions the dominant term is almost independent of H except for the In(^) 
term. Also we have neglected the F dependent terms completely due to T^ F. However, in 
a more realistic situtation the inclusion of the scattering term is necessary. 

In Fig. 121 we show Ia,b{0,4>) which determines the d,(j)- dependence of specific heat and 
(fh dependence of thermal conductivity k^z for various fixed 9. For planar field (0=90) cusps 
appear when the field points along the node directions which are equivalent to the [100] and 
[010] axis directions (0 =0, ±90°) for the A- phase. In the B- phase only twofold oscillations 
with cusps for field direction along [010] (0 = ±90°) occur. 

Very recently from the fiSH experiment a possibility for triplet sc in PrOs4Sbi2 is suggested 
If this is indeed the case, we can still reproduce the angular dependence of thermal 
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e = 45° 




-90 -45 45 90 

* [deg.] 

Fig. 3 - Comparison of angular functions Ia,b(6, (f)) for the fourfold (A) and twofold (B) gap functions. 
Inset show polar plots of the two gap functions with nodes along [100] and [010] directions for A- 
phase and along [010] for B- phase. 



conductivity as observed for both A-phase and B-phase, using a gap function |A(k)| exp(±z(p) 
which has the same modulus as eq. Q multiplied by a phase factor containing the azimuthal 
angle Lp of k. This would then be a "p+h" wave gap instead of the "s+g" wave gap in eq. 
Alternative triplet gap hmctions have been proposed by p2 23 24 . On the other hand, in this 
instance the impurity scattering is similar to d- wave superconductors. This gives a dominant 
term in k^z proportional to H2 , which is also different from experiment j^j- In principle, field 
and temperature dependence and impurity dependence of the angular part of magnetothermal 
conductivity will allow to distinguish between singlet and triplet case. For this purpose more 
detailed experiments for at lower temperatures (< 0.2 K) are highly desirable. 

Concluding remarks. - We present here a simple nodal hybrid gap function model for 
A(k) which appears to describe many features of the angular dependent magnetothermal con- 
ductivity in the multiphase superconductor PrOs4Sbi2. From the thermal conductivity we 
infer the presence of point nodes at k=(l,0,0), (0,1,0), (-1,0,0) and (0,-1,0) in A-phase and 
at k=(0,l,0) and (0,-1,0) in B-phase. We have shown that the hybrid nature of the sc order 
parameter leads to a very unusual effect of impurities. Already for arbitrary small scattering 
r an energy gap opens immediately with related low temperature exponential behaviour of 
specific heat and thermal conductivity. This implies the extreme sensitivity of s-|-g- wave 
superconductivity to the presence of impurities. We hope that the present work will stim- 
ulate further experimental investigations in the muliphase superconductivity of skutterudite 
PrOs4Sbi2. 
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